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Abstract 

The unsteady boundary layer due to a single 
rectilinear vortex filament approaching a two- 
dimensional stagnation point is investigated. 
Assuming the vortex remains far from the surface, 
incompressible potential- flow theory is used to 
determine the time -dependent inviscid flow field. 
The unsteady boundary- layer equations are solved 
by an alternating-direction-implicit finite- 
difference method. Two mechanisms which cause 
fluctuations in heat transfer are the unsteady 
velocity field in the boundary layer and secondly, 
the unsteady total temperature at the edge of the 
boundary layer. The relative importance of these 
mechanisms is dependent upon the total temperature 
fluctuations relative to the imposed temperature 
difference. As a vortex approaches a stagnation 
point it may be forced to one side of the 
stagnation line or the other, depending on its 
initial position. Results are presented for both 
of these cases. 


Introduction 


High-frequency heat- transfer fluctuations are 
a source of concern in a number of gas-turbine 
flow environments, since the thermal stresses 
which they cause can limit fatigue life even when 
the time -averaged level of heat transfer is well 
within design limits. The Space Shuttle Main 
Engine (SSME) turbopumps are one instance where 
this phenomenon has received attention, as part of 
the SSME Durability Program. The research 
reported in the present paper is part of an 
experimental/analytical effort that has been 
supported by the Durability Program, at the 
Calspan-UB Research Center (CUBRC)^ . 

In the experimental portions of the latter 
effort, time - resolved measurements of the 

distributions of heat -transfer rate have been made 
in a shock-tunnel facility, using full-scale 
turbine hardware operating under rated conditions. 
The data come from thin- film heat- transfer gages, 
which allow microsecond time- scale and millimeter 
spatial-scale resolution. 

The data taken in the experimental part of 
the effort display a very rich frequency content, 
including a number of very sharp 'spikes' which 
rise to factors as large as two above the mean 
level, and drop to values near zero, over an 
interval on the order of microseconds. The 
application of phase averaging to these data tends 
to remove most of this 'spikiness', allowing the 
fundamental blade - to-blade periodicity (with 
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periods on the order of 100 microseconds) to stand 
out. The interpretation of this fundamental 
unsteadiness in terms of laminar/ turbulent 
transition triggered by wake cutting has been 
pursued in other papers from this program 4, . 

However, the question remains as to the 
source of the higher- frequency spikes, especially 
in regard to their property of displaying heat- 
transfer levels well below the mean. The lower 
limits of these excursions are difficult to 
explain in terms of fluctuations in the thermal 
boundary -layer thickness; these might account for 
heat- transfer reductions on the order of a factor 
of two due to thickening of the thermal layer, but 
levels approaching zero exceed the bounds of this 
mechanism. 


An alternate explanation that may account for 
these down-going spikes is the passage of a 
concentrated vortex over the blade leading edge. 
These vortices, shed from the trailing edges of an 
upstream stator row, have decayed only slightly by 
the time they are intercepted by a downstream 
rotor. A number of recent studies have described 
the very rapid velocity-boundary- layer changes 
that occur during the close approach of a vortex 
to the surface®* 8 . In addition, the flow 
unsteadiness causes a momentary redistribution of 
stagnation temperature 9 . The present work was 
undertaken in an attempt to clarify whether the 
combined effects of velocity and temperature 
unsteadiness during vortex passage over a leading 
edge might explain the spikes seen in the 
experimental program. 

The net result of the research is an 
affirmative answer to this question; it does 
appear that the combined effects of the 
fluctuations in velocity and temperature induced 
by the motion of the vortex produce a heat- 
transfer 'signature' that is qualitatively similar 
to that seen in the experiments . 


The analysis that leads to this conclusion 
invokes a number of simplifying assumptions, in an 
effort to gain as much analytical understanding as 
possible. It assumes a low mach-number inviscid 
flow, makes use of incompressible potential- flow 
theory as a source of the time - dependent vortex- 
induced portions of the flow- field, and employs 
laminar, noninteracting boundary-layer theory to 
describe the viscous layer. The geometry 
considered is a two-dimensional stagnation point; 
thus the baseline flow and heat transfer that 
occur long before passage of the vortex are given 
by the classical solutions of Hiemenz 1 and 
Frossling . When the vortex comes near the 
surface, it speeds up considerably; thus the time- 


dependent contributions to the velocity and 
temperature fields become pronounced. 



The unsteady boundary - layer equations are 
solved by an alternating - direct ion- implicit 
method'* 1 , for a prescribed vortex trajectory. 
The temperature at the boundary -layer edge is 
allowed to have a constant level different from 
the wall temperature, and in addition to undergo 
the fluctuation caused by energy segregation 
during the vortex passage. At this level of 
approximation, the velocity and temperature fields 
are uncoupled, and the temperature in the boundary 
layer can be expressed as the sum of two terms, 
corresponding to the imposed temperature 
difference plus the induced fluctuation. The 
resulting heat- transfer rate is the sum of two 
corresponding terms. It is clear that very low 
heat- transfer levels with very short duration are 
possible with certain combinations of the 
amplitudes of the two temperature fields. 

Formulation 

Consider a steady two-dimensional stagnation- 
point flow; let the flow be directed toward a wall 
located at the plane y = 0, with the stagnation 
point at x - y - 0, see Figure 1. The overbar is 
used to denote dimensional quantities. The 
temperature at the wall, T w , as well _as the 
stagnation temperature far from the wall, T 0 , will 
be constant. 


The steady portion of the middle term is 
proportional to the square of the freestream Mach 
number, while the appropriate scale for the time- 
dependent potential term is the square of the 
vortex strength divided by its distance from the 
surface. Since the latter scale is not affected, 
to first order, by the freestream Mach number, it 
is legitimate to approximate the vortex motion by 
incompressible- flow potential theory. 

Because the distance of closest approach of 
the vortex is unknown in advance, there is no 
convenient length scale that can be used to define 
a dimensionless measure of the relative importance 
of time - dependent temperature fluctuations 
compared to those of steady adiabatic compression. 
The portion of the solutions below that reflect 
this unsteadiness will become large whenever the 
vortex approaches sufficiently close to the 
surface. Of course, the distance of closest 
approach must be greater than that value where the 
viscous core and the boundary layer might 
interact; the cases shown below are far removed 
from this limit, so that the model of a potential 
vortex approaching a thin boundary layer can be 
used . 

Introducing the following non-dimensional 
variables 


Now, in addition, let a single inviscid two- 
dimensional vortex with strength, V, approach from 
upstream infinity. A vortex Reynolds number may 

r 

be defined as Re v - , where v is the kinematic 
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viscosity. It is assumed that Re y is large so 
that the vortex can affect the thin viscous region 
near the wall while remaining in the effectively 
inviscid region far from the wall. 
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The motion of the vortex and unsteady 
velocities and temperature outside the boundary 
layer can be found using standard potential flow 
techniques 1 . The techniques actually used are 
those of classical incompressible flow, which may 
at first seem surprising, since they are being 
applied to a subsonic compressible - flow problem. 
The justification for this simplification is as 
follows: the temperature in the inviscid flow is 

governed by the compressible- flow energy equation: 
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where x and y are coordinates tangent and normal 
to the surface_ respectively ; also t is time, p is 
the pressure, T is the temperature and u and v are 
the tangent and normal velocities respectively. 
The constant A is the stagnation- point velocity 
gradient which determines the rate at which the 
velocity increase_s in the x direction; also the 
constants p 0 and T 0 are the steady total pressure 
and total temperature, respectively. 


The solution of this equation, for isentropic 
flow, is the compressible Bernoulli equation: 



'OD 


1 + 


i« 21 




( 2 ) 


The velocity of the vortex is 

- V» * 5^7 <‘> 
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The middle term in this equation describes the 
adiabatic-compression mechanism of heating, while 
the third term is what accounts for the time- 
dependent fluctuations in the external 
temperature. The former set tends to be spatially 
uniform in the vicinity of a stagnation point, 
whereas the latter shows significant spatial and 
temporal variations. (The middle term also 
contains a fluctuating part, in an unsteady flow). 


It can be shown that if the vortex is located 
at (x VQf y VQ ) at t-0 then the position of the 


vortex is given by 
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The unsteady tangential velocity and the 
|teady temperature induced at the wall are given 


proportional to its distance from the surface. 
The non-dimensional boundary layer equations for 
mass, x-momentum, and energy are 
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also the total temperature is given by 
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where the pressure is known from the solution of 
the inviscid problem. Also, Pr is the Prandtl 
number and E is defined by equation (12) . 


Equations (13-15) are 
boundary conditions 

at r? - 0 

as r) ° o 


subject to the following 

u - v - 0 - 0 
u - u e , 0 - $ e (16) 


The importance of the bracketed term in 
equation (10) will depend on the value of E and 
also on how large the terms within the bracket 
become during the passage of the vortex. To get 
an idea of an appropriate value for E, consider 
flow over a leading edge of radius a with 
freestream velocity u ffl . Then it can be shown that 
_ fu 1 f 

A - 0 hr 22 and — - ©(u^aA) where A is a constant 
on the a order otunity (i.e. A-0.2)®. So 
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is the Eckert number. 
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Notice this can also be written in terms of the 

u_, 

gas 


mach number M - ~ 
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where for an ideal 

c oo " J tRToo > 7 is the ratio of specific heats, R is 
the gas constant. Then 
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To assume incompressible flow, as has been done in 
the present formulation, the double limit as M 2 -> 


0 and 



0 must be taken. 


From the above equation one can see that it is not 
inconsistent to choose a finite value of E. 


where u and 
respectively . ' 


are given by equations (8) and (9) 


To facilitate the numerical solution of 
equations (13) and (14) an unsteady stream 
function is defined by 
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It may be shown that equations (17) satisfy the 
continuity equation identically and the x-momentum 
equation can be written 
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The freestream velocity is often included in 
the definition of the stream function to simplify 
the boundary condition. In the present case this 
is not convenient since the unsteadiness of u 

e 

will cause it to vanish at x-positions where the 
velocity profile is not identically zero. 

Owing to the linearity of the energy equation 
(15) , it is possible to define two new temperature 
fields such that 

$ - + E$ 2 ( 20 ) 


Since both A and E c are expected to have 
values somewhat less than one, values of E which 
are quite small are investigated. Results are 
presented for E - 0.0 and E - 0.02. 

Even for a small value of E, the bracketed 
term in equation (10) can be important during the 
passage of the vortex. The major effect of the 
vortex occurs directly beneath the vortex where 
(x-x v ) is small, so in general the fluctuation in 
total temperature is directly proportional to the 
tangential velocity of the vortex and inversely 
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subject to the boundary conditions 
at r? - 0 0 1 “ 0 2 “ 0 


with boundary conditions 

at r)-0 " 0 


as r) -*■ oo -► 1 , -+ (23) 
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Notice that a parametric investigation of the 
effect of E on 0 requires solving for and d 2 
only once. 

Initial Condition 

To begin the numerical calculation the 
solution over the entire (x,y) plane must be known 
at the initial time t-o. It is possible to find 
an asymptotic solution, valid when the vortex is 
far from the wall, which can be used for the 
initial condition, as follows: 

Expanding in powers of (l/y v ) the inviscid 
results for u & and d & become 
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By expanding the unsteady stream function and 
temperatures as 
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and introducing these relations into equations 
(18) , (21) and (22) the following set of ordinary 

differential equations results 
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Note that equations (31) and (33) with the 
boundary conditions of equation (37) correspond to 
the steady stagnation-point problems of Hiemenz 1 
and Frossling . Each of equations (31) *(36) 
together with boundary conditions (37) are solved 
numerically using methods which are second— order 
accurate . 

Assuming the initial vortex position is 
chosen large enough, equations (28) -(30) will 
produce an accurate initial condition for the 
problem. 

Numerical Method 

The numerical method used to solve the 
boundary-layer problem described by equations 
(18), (21) and (22) is similar to that described 

by Peridier and Walker (1988)'* . The method is 

a factored alternating-difference- implicit method 
(ADI) which uses upwind -downwind differencing for 
the first-order derivatives. The following is a 
brief description of the method. For a more 
complete description see Peridier and Walker 
(1988) 7 * 12 . 

A rectangular grid is defined, where the 
indices i,j,n denote the x-direction, the re- 
direction and time, respectively. Constant grid 
spacing is used, where Ax, A r? and At denote the 
respective step sizes. Consider in general the 
following equation 
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1 

shows 

how 


represent each of equations (18), (21) and (22). 
A finite-difference approximation to equation (38) 
is written at the point (i,j,n+l/2). The 
coefficients are approximated using a simple 
average, for example 
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where the overbar here denotes that averaging has 
been done. The derivatives are also averaged in 
this manner. The second derivative is 
approximated by a central difference. The first 
derivative terms are approximated using an 
upwind-downwind scheme according to the sign of 
the coefficient. For example 
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Equation 40 may be written more compactly by 
defining the operators and as 

X+ _ 6 gSgn(Pj j )/2 x - _ 6 E - S gn(P t -j)/2 (41) 
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where 8 ^ is the central difference operator 
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from the previous iterate. The procedure is 
repeated until the maximum difference between 
iterates is below a given tolerance. Notice that 
for a linear equation, as is the case for the 
energy equations (21), (22), the coefficients need 
not be iterated upon. However, as will be 
discussed below, the application of the downstream 
boundary condition may require iteration. To 
update the coefficients in the case of the 
momentum equation (18) , the tangential velocity u 
must be integrated to get xj>, This is done using 
Simpson's 1/3 rule 1 . 


and is defined by 


Boundary Conditions 
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Boundary conditions are applied around the 
(43) entire perimeter of the numerical domain. At the 

surface, i.e. , r?-0 


It follows that equation (40) can be written 
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As y) -*■ «> the solution must asymptotically approach 
the freestream values. The freestream boundary 
conditions are applied at r) - *? max where y max is 
chosen large enough so that increasing y max does 
not affect the solution. A typical value of y max 
is ~ 9. 

The x-domain is Hw x^] where is 

chosen large enough so that increasing x^^ does 
not affect the solution. The effect of the vortex 
varies inversely with the distance from the vortex 
so the boundaries of the x-domain are not expected 
to experience much unsteadiness. This fact makes 
it possible to use a first-order-accurate 
extrapolation at the x-boundary without 
introducing any appreciable error in the region 
near the vortex. The method is developed by 
writing a Taylor series expansion about the grid 
point just inside the boundary. For example, if 
x ( 1 ) - - x(2) - - Xjj^ + Ax and so on, then 


The finite difference approximation to equation 
(38) , written in factored form is 
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This method is second order accurate in space 
and time. The solution procedure is as follows. 
Equation (47) is solved along all lines of 
constant r? . Then equation (48) is solved along 
all lines of constant x. The solution of equation 
(47), as well as equation (48), is accomplished 
using a direct method of solution for tridiagonal 
problems. The method involves a combination of 
cyclic reduction 1 and the Thomas algorithm 1 . If 
equation (38) is nonlinear several iterations may 
be required to advance the solution in time. The 
coefficients are first estimated using previous 
time— plane values and then updated using values 
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using values from the previous iteration. 
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where the prime ( ) ' denotes the previous 

iterates. Thus the boundary condition at x - x max 
is 


F(1 , j ) - F(2 , j ) - - |[(F'(3,j) - F'(l.j)] + 0 (Ax 2 ) 
L J (52) 

Similarly at x - x max 
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+ 0(Ax 2 ) (53) 

Notice that by applying the boundary 
conditions in this way more iterations may be 
required. Typically only 2 iterations are 
required for the endpoints to come within the 
tolerance . 


5 


TABLE 1 


Equation 

F 

T 

P 

Q 

M 

18 

u 

i 

di> 

dx 

— u 

_ ^ 
3x 

21 


1/Pr 

—v 

-u 

0 

22 

*8 

1/Pr 

-v 

-u 

fdul 2 3p dp 

— + — + u — 

3t dx 


Results 

It is possible to investigate all cases by 
choosing the sense of the vortex to be 
counterclockwise, since changing the sense of the 
vortex and the sign of x VQ would result in a 
mirror image of the previous results. Two 
representative cases are presented. Each case is 
determined by the initial vortex position 
(x , y ) . For each of the two cases heat 
transfer results are presented for E-0 . 0 and 
E-0.02. The Prandtl number is held fixed and 
taken to be Pr -0.72. Notice that since the 
vortex strength, T, and the stagnation point 
gradient, A, are included in the 
nondimensionalization they do not appear 
explicitly in the calculation. 

For the first case the initial vortex 
position is (x v , yy 0 ) " (0.025,20.0). This is 
referred to as the right case, since as the vortex 
approaches the surface it is pushed to the right 
of the geometric stagnation point, refer to Figure 
2(a). In Figure 2(a) the vortex path is shown and 
three points are marked corresponding to t - 1.50, 
3.00 and 3.40 as shown. Results are then 
presented as functions of x at these three times. 
As noted earlier, the sense of the vortex is taken 
to be counterclockwise. The vortex causes an 
increase in the external velocity, Figure 2(b), 
and its effect is felt mainly beneath the vortex. 
The resulting wall shear stress is shown in Figure 
2(c). The wall shear stress increases beneath the 
vortex, but as the vortex comes nearer to the 
surface, the shear stress takes on values below 
the steady-state value in the region in front of 
the vortex (i.e. x > x^,) . This occurs because the 
pressure gradient in this region is not as 
favorable as the steady value. At later times the 
pressure gradient would become adverse in front of 
the vortex and cause flow reversal, eventually 
causing a breakdown of the boundary-layer 
approximation. Figure 2(d) shows the displacement 
thickness as defined by 



JA 


There is a tendency for the boundary layer to thin 
behind the vortex and thicken in front of it. 
Note that near x-0 the displacement thickness 
fails to give meaningful information, since it is 
undefined when u g — 0. Plots of the 
streamfunction , Figure 3, also give the indication 
of thinning behind and thickening in front of the 
vortex. 


Figure 4 shows the external total 
temperature, the thermal boundary layer thickness, 
and the resulting heat transfer for E-0.0 and for 
E-0.02. The thermal boundary layer thickness is 
defined by 


* x - — L “ «? 
u 

\ A 

Figure 4(a) shows that the external total 
temperature for E-0.02 increases by as much as 30% 
below the vortex. For E-0.0 the external total 
temperature is identically equal to one. In 
Figure 4(b) 6 T is shown only for E-0.02 because 
the differences between E-0.00 and E-0.02 are very 
slight. Note that as with the displacement 
thickness the trend is to thin behind the vortex 
and thicken in front of it. 

Figure 4(c) shows heat transfer results for 
E-0.0 and E-0.02. The dimensionless heat transfer 
is defined by 


0-. 990 . 


(55) 


where 


39 } 


3 1?. 


l»l-0 


-k(T 0 -T w )l 



V 

> I 

t 


- 3T 
- k — 

*y. 


y-0 


(56) 


(57) 


Notice that for E-0.0 the heat transfer is 
higher where the boundary layer is thinner and 
vice versa. For E-0.02 there is an additional 
increase in heat transfer due to the increased 
external total temperature. Figure 5 shows the 
heat transfer as a function of time for x - -2 , 
-1,0, 1,2. From this it can be seen that the 
vortex passage has little effect on negative 
values of x and even at x-0 . It is interesting to 
note that a time average of the heat transfer 
would result in a net increase for this case. 
This is most easily seen at x-2.0. Since the time 
is non- dimens ionalized by A, and A may take on 
very large values , the time scale in Figure 5 is 
greatly exaggerated. 


The second case corresponds to an initial 
vortex position of (x^.y^) - (-.185,20.0) and 
will be referred to as the left case. Figure 6(a) 
shows the vortex path for this case. As in the 
right case results are presented as functions of x 
for the three times corresponding to the points 
marked in Figure 6(a). For this case, the vortex 
is pushed toward the left (i.e. x<0) where the 
steady external velocity is negative, Figure 9(b). 
As in the right case the vortex causes a positive 
addition to the external velocity which drives it 
towards zero . When the vortex gets close enough a 
region of " inviscid flow reversal" occurs beneath 
the vortex. As can be seen in the plots of wall 
shear, Figure 6(c), and streamfunction, Figure 7, 
a region of flow reversal forms beneath the 
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vortex. This however, does not imply a breakdown 
of the boundary layer approximation since the flow 
is unsteady Note that the displacement 
thickness becomes singular where the outer 
velocity goes through zero and thus loses its 
usefulness for the present case . As a qualitative 
check on the boundary layer approximation, the 
maximum normal pressure gradient was calculated 
numerically from the results. This is done by 
using the y-momentum equation which was neglected 
in arriving at the result: 


dp dv dv dv d 2 v 

— - — + u — + v 

dr? dt dx dr? dy 2 


(58) 


The right case which showed no signs of separation 
was found to have maximum normal pressure 
gradients comparable to the left case . Flow 
separation would be expected to occur eventually 
for either case. 


Figure 8 shows the external total 
temperature, thermal boundary-layer thickness, and 
heat-transfer results for the left case. Note 
that the external total temperature (Figure 8a) is 
decreased as the vortex passes. The thermal 
boundary layer thickness (Figure 8b) , behaves 
similar ,to the right case, i.e. thinning for x < 
x v and thickening for x > x v> Figure 8c shows the 
resulting heat transfer. As expected it is higher 
where the thermal boundary layer has thinned and 
lower where the thermal boundary layer is thicker 
for E-0.0. For E-0.02 however, the heat transfer 
is lower than the E=0.0 result because of the drop 
in external total temperature . The time traces of 
heat transfer, Figure 9, show that values of x > 0 
and even x - 0 are only slightly affected by the 
passage of the vortex. For this case there 
appears to be a net decrease in the heat transfer. 

For each case the numerical domain was 0 < rj 
<9 , -6 < x < 6 , with constant spatial grid sizes 
of A r? - .125 and Ax = .05. The time step was 
At “ .01 for t < 2.5 and .005 for t > 2.5. Each 
case took approximately 300 seconds of CPU time on 
a CRAY-XMP . 


Concluding Remarks 


departures from two - dimens ional i ty , to name just a 
few. Thus, it is not possible to make a direct 
comparison with the measurements. However, 
further study of these effects can now be carried 
out with confidence that the basic mechanism is a 
significant one in producing heat - transfer 
fluctuations in the stagnation region. 
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